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Abstract 

Basic physical principles are considered that are responsible for the 
origin of dynamic air flow upon condensation of water vapor, the partial 
pressure of which represents a store of potential energy in the atmosphere 
of Earth. Quantitative characteristics of sucli flow arc presented for several 
spatial scales. It is shown that ma^cimum condensation-induced velocities 
reach 160 m s~^ and are realized in compact circulation patterns like 
tornadoes. 

1 Introduction 

Atmospheric air on Earth conforms to the ideal gas law with a high 
accuracy. The main physical property of the ideal gas is that its 
equation of state does not depend on molar masses of the gas mix- 
ture constituents. At fixed temperature, a given value of pressure 
can be obtained for a mixture with arbitrary molar masses by set- 
ting the value of molar density of the gas. The second important 
property of atmospheric air is the presence of a constituent that un- 
dergoes condensation under terrestrial temperatures and pressures 
- water vapor. Molar density of the moist air mixture is equal to 
the sum of molar densities of the dry air and water vapor. 

In a motionless atmosphere (with the impact of the greenhouse 
substances neglected) the air temperature would be the same at 



all heights. Vertical distributions of all gases would follow the hy- 
drostatic Boltzmann's distribution according to their molar masses. 
The exponential scale height of water vapor with molar mass 18 g mol 
would be about 13 km, while the scale height of the major air 
constituents, nitrogen and oxygen, with similar molar masses ~ 
30 g mol~^ would have scale heights of approximately 8 km. In 
this case water vapor would be saturated at the surface only (due 
to the contact with the liquid hydrosphere) and had an undersat- 
urated concentration elsewhere. There would be no evaporation or 
condensation in such an atmosphere. 

However, irrespective of the presence or absence of the greenhouse 
substances, such a hypothetical distribution of vertically isothermal 
moist air appears to be unstable. Any fluctuation leading to an 
upward displacement of an air volume results in adiabatic cooling 
of the rising air. Air temperature drops such that the equilibrium 
water vapor concentration dictated by Boltzmann's distribution be- 
comes oversaturated at all heights where the air ascends. This causes 
the water vapor to condense. Its concentration decreases down to 
the saturated concentration. Condensation diminishes the total air 
pressure and disturbs Boltzmann's distribution of moist air. The 
vertical gradient of air pressure becomes greater than the weight of 
a unit air volume. There appears an upward-directed force acting on 
a unit air volume. Static equilibrium of moist air in the gravitational 
fleld is no longer possible. There appears a rising flow of air masses 
induced by condensation. The process of condensation is sustained 
by continuous evaporation of water vapor from the hydrosphere. 
The upward-directed force that acts on moist air causing it to rise 



adiabatically was termed the evaporative-condensational force ( 


Gor- 


shkov and Makarieva, 


2006 


Makarieva and Gorshkov 


2007, 


2009a 


)■ 



Most part of the condensed water vapor leaves the atmosphere 
via precipitation. A minor part is maintained in the atmosphere by 
the rising air flow. This imposes a drag force on the rising air flow 
and leads to a reduction of the vertical velocity. However, unlike 
the flow of water vapor which condenses as it rises, the flow of the 
dry air components which conserved their mass, cannot be unidi- 
mensional (vertical). There inevitably appear horizontal legs in the 
condensation-induced air circulation. Condensation of water vapor 
in the ascending air produces both vertical and horizontal pressure 
gradients. Therefore, the presence of water vapor in the atmosphere 
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contacting with a liquid hydrosphere leads to the formation of three- 
dimensional circulation patterns. 

In the following sections we derive relations between the verti- 
cal and horizontal components of air velocities and condensation- 
induced pressure gradients. We further use the obtained results to 
describe large-scale circulation with approximately constant veloci- 
ties when the pressure gradient force and the turbulent friction force 
coincide. We also apply these relations to describe hurricanes and 
tornadoes, where the pressure gradient forces appear to significantly 
exceed the turbulent friction forces. So far the condensation-induced 
air motions have not received a consideration in meteorology. Sev- 
eral relevant observations regarding the conventional approaches are 
made in the footnotes to the main text. 



2 Continuity equation for moist air 

The equations of state for moist air as a whole, as well as for its 
components - dry air and water vapor, include one and the same 
universal molar gas constant R and do not depend on molar masses 
and mass densities of the component^ 

p = NRT, p, = N,RT, pd = NdRT, (1) 

where p, N, p^, Ny, p^, are the pressure and molar density of 
moist whole, water vapor and dry air, respectively. 

In a circulating atmosphere where no condensation takes place, 
the ratios of molar densities of all components at all heights affected 
by the are equal to their mean atmospheric values. The reason is 
that the diffusional velocities that would restore Boltzmann's distri- 
butions depending on molar densities and molar masses of the com- 
ponents, are small compared to the dynamic velocities of the air flow. 
According to observations, mixing ratios of the non-condensable air 
constituents and the molar mass of dry air are the same at all heights 
in the troposphere. 

In the absence of condensation the ratio 7 = N^/N should not 
change at any changes of pressure and temperature. Consequently, 

^Note that this fundamental universality of ideal gas is masked in the meteorological lit- 
erature by the common usage of mass densities p = MN, = MyN and p^ = M^N along 
with mass gass constants Rair = R/M, R^ = R/My and R^ = R/Mj. Usually Rair ~ Rd is 
denoted as R, while the universal molar gas constant is practically never used (e.g. , [Glickman[ 
|2000[ |. 
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the process of condensation should be reflected in the changes of ra- 
tio 7. The molar rate of condensation per unit volume (mole m~^ s~^) 
that is caused by the decrease of air temperature with height z due 
to the adiabatic ascent of moist air with vertical velocity w, is equal 
to 

^^1^ = ^M^-T7^ ' (2) 



dz \ dz N dz ^ 
_ Njj d'y dN^ dN dp^ dp ,^ 

The meaning of Eq. ^ is physically transparent. Condensation 
rate is determined by the change of water vapor concentration minus 
the change of total air concentration that is not related to conden- 
sation. We emphasize that the second term in brackets in Eq. ^ 
comprises the relative change of molar density of moist air as 
a whole rather than molar density Nd of its dry component. This 
reflects the fact that restoration of equilibrium pressure distribution 
upon condensation affects the air mixture as a whole, including the 
remaining water vapor. 

Since the condensation rate ([2]) is a function of molar (not mass) 
densities, it is convenient and physically transparent to write the 
continuity equation for ideal gas in terms of molar densities as well. 
In Cartesian coordinates (assuming that there is no dependence of 
the flow on y) the continuity equation takes the form, see Eq. ([2]): 

dNu dNw fdN, N^dN\ 

+ ^^ = ^\^-^^]- (4) 



dx dz \ dz N dz 

Taking into account that N = N^ + N^ and 

dNdU dNdW 



dx dz 

we have from Eq. Q: 



0, (5) 



dN^u dN^w fdN^ N^dN\ 



dx dz \ dz N dz 

Expanding the derivatives in Eq. ([6]) and multiplying both parts of 
the equation by N/N^ we obtain 

dN^ ^dN ^fdu dw\ _ ^ 
N^ dx dz \dx dz ) ' 
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Now expanding the derivatives in the left hand part of Eq. ^ and 
using Eq. ([T]) we obtain 

""[dx Ny dx J '"^[dz N dzj- ^ ^ 

Using definition ^ we can re-write the continuity equation ([s]) as 

^91n7 ^ /gN 

dx dz 

Now, turning from molar densities N and to pressures p and py 
([T]) we can see that temperature dependencies cancel from Eq. ^ 
and the latter become^ 

=wp'y (——--— ^ (10) 
\dx dx J \py dz pdzj 

According to Clausius-Clapeyron equation, saturated pressure py of 
water vapor depends on temperature only. When the considered 
area is horizontally isothermal (temperature T does not depend on 
x) we have dpy/dx = 0. This condition presumes the existence of an 
inflow of water vapor caused by local evaporation from the hydro- 
sphere. This inflow partly compensates condensation that occurs in 
the atmospheric columi]|^ In this case Eq. (10) takes the form 

dp w Ap 
dx u /i^ 



Ap = p7 = p„ h-' = h-'-h-\ (11) 



"7 



1^ = 0. (12) 

Py dz p dz dx 



■^Equation (j4j is equivalent to the sum of Eqs. |5| and |6|. If there is no condensation, the 
right-hand part of Eq. |4]l is zero. Equations |4jl, ([sj^and |6pcan then be re- written in terms of 
mass densities p = N/M, = N^/M^i and p„ = N/M^. Here My and Mjj are constant (M^j is 
constant in agreement with observations because the molar ratios of the dry air constituents 
do not change with height). But the molar mass of moist air M = ■yMy -|- (1 — 7)M^ = 
^d{^ ~ 0.387) depends on 7 and, consequently, on z. Therefore, if the condensation rate in 
the right-hand part of Eq. Q is not zero, it will change substantially upon transition from 
molar to mass densities. Specifically, if condensation rate is written in form of the right 
hand part of Eq. ^ with N and Ny changed to p and pv , respectively, this will lead to the 
appearance of an incorrect multiplier M^/M in the right hand part of Eq. This would 

contradict the physical meaning of the ideal gas equations of state (111. 

^It is assumed in the conventional meteorology that the cause of atmospheric circulation 
is the fact that the surface is not horizontally isothermal due to external differential heating. 
The condensation-induced circulation, in contrast, arises on a horizontally isothermal surface 
(although this is not an indispensable condition). Horizontal temperature inhomogeneities 
that can be observed after the circulation has established are the consequences of the horizontal 
inhomogeneity of the process of water vapor condensation. 
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All magnitudes entering (11) and (12) depend on x and z. Height 
has the meaning of characteristic height where all water vapor 

condenses. Heights and h are the scale heights of water vapor 

and moist air, respectively. 

In a large-scale stationary circulation, where constant friction 



forces compensate the equally constant pressure gradient forces (Makarieva 



and Gorshkov, 2009a), horizontal velocity u does not change with 



X. Velocity w should be understood as the vertical velocity av- 
eraged over height h^. The flux of air enters the circulation area 
horizontally via a vertical cross-section of area Dh^ and leaves the 
circulation area vertically across a horizontal cross-section of area 
DL. Here L is the horizontal dimension (length) of the circulation 
area counted along the x-axis, D is the circulation width counted 
along the y-axis perpendicular to the horizontal air flow. Taking into 
account that the number of air mols riin that enter the circulation 
area differ from the number Uout of air mols leaving the circulation 
area by a relatively small number of mols of condensed water vapor, 
— Tiout ~ 7^m ^ ^in, to the accuracy of 7 <^ 1 we can write (see 
also footnote 6 below): 



Dh^u 



DLw, 



w 
u 



L ' 



Putting (13) into (11) we obtain 



dp Ap 
dx L 



, Ap = p-f = p^. 



(13) 



(14) 



Thus, total horizontal pressure drop is equal to partial pressure of 
water vapor. 

To arrive to Eq. ( [I4| ) three equations have been used: Eq. ^ for 
Ny, equation dp^/dx = 0, see (12), that reflects that the circulation 



area is horizontally isothermal and that water vapor is saturated, 
and Eq. ([s]) for molar density N^. The latter equation does not 
include the condensation rate term (|2]), which is present in ^ and 
([6]). Equation (13) namely arises from Eq. ([s]) for A^^^ ^ N^. 

It follows from Eq. (14) that p^ = Ap represents a store of po- 



tential energy, which is converted to the kinetic energy of moving 
air masses as the water vapor condensed By integrating (14) we 



* In the absence of friction at constant pressure gradient the horizontal velocity grows with 
distance x as prescribed by Bernoulli's equation; = Ap = 7p. Due to 7 <g; 1 we observe 
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Figure 1: Air streamlines for a large-scale circulation of length L and height hj. 
The continuity equation ( 13 ) is uh^ = wL. 



obtain the following expression for potential energy P{x) = p{x) 
that is defined to the accuracy of a constant term: 

Av 

P{x)=p{x)=p{0) + -f-x. (15) 



In the incoming air flow at x = L water vapor is present everywhere 
in that part of the atmospheric column that is affected by the cir- 
culation. At X = this water vapor has been completely used up 
and is partially replaced by locally evaporated water vapor. 

Strictly speaking, Eqs. (13) and (14) are valid for a bunch of 



streamlines that enter the circulation area in the horizontal direction 
and leave it in the vertical direction. The number of horizontal 
streamlines present in the atmospheric column at a given x decreases 
as one travels inside the area at the expense of those streamlines 
that have left the area in the vertical direction. Velocitites u and 



w (13) and pressure gradient (14) are the same for each streamline; 



that Ap is much smaller than p. On the other hand, we have Ap/p = 7 = Ap/p- Mass density 
p changes little over the distance where air pressure changes by Ap: Ap = ■yp <^ p. From this 
in the meteorological l iterature it is concluded that changes in density p due to condensation 
can be neglected (e.g., [Sabato[[2008[ l. . The resulting physical inconsistency is masked by the 
fact that changes of p and the formation of Ap are conventionally ascribed to independent 
physical causes. In particular, it is common in circulation models to take pressure profiles 
from observations. The relative smallness of horizontal pressure drop Ap is never discussed. 

But if we put Ap/p = Ap/p equal to zero, no velocity can form and the circulation cannot 
exist. This statement is general and does not depend on why Ap and Ap = po — p actually 
form. Mathematically, the error can be spotted as follows. From the equation of state we 
have p = Cp, where C = RT/M = const, for the considered horizontally isothermal surface. 
From Bernoulli's equation we then have = 2Ap/p = 2C'Ap/p = 2C{Ap/po){l+Ap/po + ...). 
As one can see, discarding Ap compared to p does indeed correspond to discarding the term 
of a higher order of smallness. But with respect to the pressure gradient, the main effect is 
proportional Ap, which is the term of the first order of smallness. If the smallness is set to 
zero, the effect disappears. 
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the mean horizontal velocity in the column linearly decreases, while 
the height of ascent linearly grows, Fig. 1. All the above forms of 
the continuity equation, Eqs. (11)-(14), describe the kinematics of 
the air flow at given values of u and w. These velocities should 
be determined from Euler's equations with an account of frictiorj^ 
Thus, condensation occurs as the moist air ascends in the vertical 
direction and is maintained due to evaporation from the horizontal 
Earth's surface. The associated condensational-evaporative force 
induces makes the moist air masses circulate along the streamlines 
that include both vertical and horizontal regions. 



3 Condensation in the adiabatically ascending 
air 

According to Clausius-Clapeyron equation, we have 



Pv T 

1 dpy r dT 

-^- = h^ , r = - — , ^ = —, ^0 = 



(16) 



where = 45 kJ mol ^ is the molar heat of vaporization (latent 
heat). 



^The major friction force per unit area of the Earth's surface, which opposes the pres- 
sure gradient force that would otherwise accelerate the air, is the friction force that can be 
called gravitational, as it is proportional to the weight of atmospheric column fipgh = fips, 
fi = zx/h ~ 10~^. Here pg is surface pressure, zx ~ 0.1 m is the surface roughness (it is 
proportional to the height of vegetation cover, oceanic waves etc.), h ^ 10 km is the scale 
height of the atmosphere. The same form of gravitational friction, ppgh, is due to frictional 
dissipation of liquid drops precipitating or suspended in the atmosphere. In this case p has 
the meaning of the relative volume occupied by the drops, p ~ ■y{w/wi,) ^ 10""", where Wf, is 
the mean downward velocity of the drops, w/wt ~ 10"'^. The gravitational friction force does 
not depend on velocities u or w and can be represented as ppgh = pUg^,, where p = Ug^/gh 
is Froude's number, tig, has the meaning of rotation velocity for the turbulent eddies that 
bud from the main air flow due to gravitational friction. The gravitational friction force pUg^ 
exceeds by 30 times the force of aerodynamic friction C£)pv? = pu^ ("g* ~ SOu^) that is 
usually taken into account in the Navier-Stokes equations and the formulation of Reynolds 



stress ( [Makarieva and Gorshkov] |2009"a[ l . 
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Moist air obeys the hydrostatic equihbrium distributio]:]0 

Functions T{z), T{z) and 7(z) can be found from the first law of 
thermodynamics under condition that the ascent is adiabatic, using 
the definition of 7 (|3| and the Clausius-Clapeyron equation (16): 



dz N dz ^ dz ' J dz dz pdz ^ 

(18) 



The first equation in ( 18 ) is the first law of thermodynamics for an 
adiabatic process, with the third term describing condensation, see 
([2]); Cp = + R = {7/2)R is the molar heat capacity at constant 
pressure. 



System of equations (16)-(18) can be re- written in a closed form: 



-f.r(.)^r.^(i- 0.387), (19) 
1^7 , r(^) M,(i- 0.387)^ . 

-fdz ^ ^ T RT ■ ^ ^ 

T RT R 

= fi-r- = 9.8Kkm~\ ha = -r-;—, /i = — . 
hd Mag Cp 

Here is the value of T{z) at 7 = 0, i.e. it is the adiabatic lapse 
rate of air temperature for dry aii[^ Eq. (19) at 7 ^ 1 coincides 



with the well-known expression for moist adiabatic lapse rate of air 



^Equation \l7\ should be more appropriately referred to as the equation of aerodynamic 
equilibrium in the gravitational field of Earth. Indeed, this equilibrium does not correspond 
to a static Boltzmann's distribution of gases with different molar masses. It arises in the 
result of the air ascent which occurs with a sufficiently high velocity w, which is the same 
for all gases (including the remaining water vapor) despite their different molar masses. In 
the system of rest of the ascending air, Eq. ^17\ corresponds to a hydrostatic distribution of 
air with a constant molar mass M and departs from the latter only insignificantly due to the 
fact that 7 <C 1 decreases with height. Vertical turbulent flux of water vapor associated with 
evaporation from the surface leads to the fact that the vertical velocity of water vapor Wv 
in the ascending air flow is always larger than the vertical velocity of air w, Wy = w + Aui^, 
| |Makarieva and Gorshkov| |2007^ . At 7 <^ 1 this does not change the flux of moist air as a 
whole, Nw = Niiw + N^Wy = Nw(l + •yAwy/w) ^ Nw. The ratio between vertical velocities 
of water vapor and air as a whole reaches its maximum value Aw^/w ^ 1 in a large-scale 
circulation where w is small. In compact intense circulation where condensation rate greatly 
exceeds the local rate of evaporation, we have Au!„ <g; w and Wy 

^In the meteorological literature one can find statements (e.g., [P5schl[ [2009[ p. SI2436) 
that the release of latent heat Ly upon condensation warms the air, as well as that latent 
heat is a major source of energy for some types of circulation (e.g., see discussion and refer- 
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temperature (Glickman, 2000) Solutions of Eqs. (19), (20) are 
shown in Fig. 2. 

Due to the large value of the ratio L^/R = ^ 5300 K, the 
dimensionless ratio ^ = T^/T is always much larger than unity. 
Therefore, at 7 — )■ 1 (this would happen if at constant temperature 
the dry component were largely removed from the atmosphere and 
the water vapor partial pressure became the dominant contributor to 



total air pressure) the adiabatic lapse rate (19) ceases to depend on z 
and tends to T -> (rd)(0.62//iO = 1-2 K km-\ Scale height K ^ 
of saturated water vapor tends to the hydrostatic equilibrium value 
K = T/{TO ^ /id/0.62 = K,tat = RT/{M,g) = 13.5 km. The 



right-hand part of the last equality in Eq. ( 18 ) turns to zero; 7 ceases 



to depend on z. The third term containing d^/dz in the equation 



of the second law of thermodynamics (18) vanishes. Condensation, 
evaporation and the condensation-induced air circulation all stop. 

The transition to hydrostatic equilibrium during adiabatic ascent 
of moist air at 7 — )■ 1 implies that the weight of saturated water va- 
por column of unit area at any height becomes equal to saturated 
vapor pressure at this height: the evaporative-condensational force 
disappears. Consequently, the adiabatic ascent of moist air can no 
longer be maintained by this force and cannot arise spontaneously. 
In this case, in contrast to the case of 7 ^ 1, the stable stationary 
state of the atmosphere is not the state of self-sustained air circula- 
tion, but the the state of hydrostatic equilibrium of motionless air 
(water vapor) with constant air temperature at all heights. 

The condensation-induced circulation of the modern atmosphere 
occurs at a global mean surface temperature that is optimal for life. 
This temperature fixes the saturated concentration of water vapor in 



accordance with the Clausius-Clapeyron equation (16). Therefore, 
the condensation-induced circulation only becomes possible due to 
the high concentration of atmospheric nitrogen, which can be called 
an air "ballast". Indeed, unlike CO2, nitrogen is not a greenhouse 



ences in |Makarieva et al.[ |2010a[ |. In reality, according to the Clausius-Clapeyron equation, 
condensation occurs when there is an external cause of cooling (e.g., adiabatic ascent). Con- 
densation cannot warm the air to a temperature higher than the air had prior to condensa- 
tion. During condensaton, the water vapor concentration decreases, which, according to the 
Clausius-Clapeyron law, correspond s to a dro p (not a rise) of air temperature. 

* We emphasize that 7 in Eqs. \20\ is defined as in Eq. |3t. At ~^ we have 

7 — >■ 1 and d^/dz — > 0. In the meteorological literature formulae l |19[ l and ( |18[ l are written 
for 7^ = Pv/Pd instead of for 7 = Pv/p |3jl. At 7 <S 1 we have 7 7£j, but at — > we have 
7£j — > 00, which leads to a physically meaningless expression d'-^^./d^ — > 00 if one replaces 7 
by 7d in the continuity equation (|9| and Eqs. ( |20| and (|2|. 
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Figure 2: Dependence on height z of the major distribution functions of moist 
air for several different values of surface temperature Tg satisfying the system 



of equations (19)-(20). 



1 Up 



(a): h(z) ^ = — ^ is the scale height il7h of the moist air pressure; 
poz 



1 

(b): h^{z)^^ = — ^ is the scale height (|16|) of the partial pressure of 

Pv oz 

saturated water vapor; 



1 d-f 



(c) : — — — is the scale height of water vapor condensation ( 20 1 ; p and py 

J oz 



are pressures of moist air and saturated water vapor, respectively. 

(d) : 7(2:) = Py{z) / p{z) is the mixing ratio of saturated water vapor; 

(e) : r(z) = —dT/dz is the moist adiabatic lapse rate of air temperature, the 
dashed line indicates the dry adiabatic lapse rate Td — 9.8 K km^^; 

(f) : ^{z)/hj{z) = —dj/dz is the relative intensity of water vapor condensation. 
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gas; unlike oxygen, gaseous nitrogen is chemically inert and cannot 
lead to excessive oxidation or fires in the biosphere. The major role 
of atmospheric nitrogen consists in making the value of 7 small at 
a fixed temperature. The condensation-induced circulation makes 



the biotic pump of atmospheric moisture (Makarieva and Gorshkov 



2007) and a hydrological cycle on land possible, thus ensuring that 
land is habitable for life. This suggests that the existing concentra- 
tion of atmospheric nitrogen should have been formed by the time 
while life started colonizing the land and the continental forest cover 
developed. 

Note that if 7 — )■ 1 due to increasing temperature (and not due to 
the removal of the dry air component), then the exponential growth 



of py = 7p in accordance with the Clausius-Clapeyron equation (16) 
makes the evaporative-condensational force grow exponentially as 
well, such that the condensation-induced atmospheric circulation is 
ensured at any small share of the dry component in the total air 



pressure, see (12), (14) and Fig. 2f. 

For a given value of 7(0) = 7s ^ 1 at the surface, 7(2;) declines 
rapidly with increasing z, Fig. 2d. Accordingly, the adiabatic tem- 
perature lapse rate T(z) ( [l9) ) rises from its minimal value at the 
surface approaching the dry adiabatic lapse rate at large heights. 
Fig. 2e. Any circulation pattern includes areas where the air masses 
rise (in the region where evaporation and condensation are more in- 
tense) and areas where the air masses descend (in the region where 
evaporation and condensation are less intense). In the region of 
descent the air masses warm adiabatically while descending, so con- 
densation cannot occur. Consequently, neglecting horizontal mixing, 
the adiabatic temperature lapse rate in the region of descent should 
be equal to the dry adiabatic lapse rate. Air temperature T and air 
pressure p related to temperature by the equation of state ([T| de- 
crease with height more slowly in the region of ascent than they do 
in the region of descent. Pressure difference Ap between the regions 
of ascent and descent has a negative value at the surface, then de- 
creases by absolute magnitude with growing height, approaches zero 
and changes its sign at a certain height Zc ~ h^, where horizontal 



velocity u changes its direction (Makarieva et al. , 2010b). 

The equation of state ([!]) can be written in the following form 

p = pgh, p = NM, h = RT/{Mg). (21) 
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This relationship for p does not in reahty depend on either M or 



which actually cancel in (21). But written in this form, presure has 
a simple physical meaning of potential energy in the gravity field, 
where h represents the height of atmospheric column, pg represents 
weight of a unit air volume in the air column, while p represents 



the total weight of the air column of a unit area. Relationshp (21) 
holds for any height with p and T depending on z. At the Earth's 
surface z = 0; surface values of all variables are denoted by low index 
s. The difference in surface weights between two hypothetical static 
air columns - in one of which saturated water vapor has condensed 
at all heights, while in the other no condensation took place - is 
equal to the difference of the right-hand parts of the first equality 



in (21) written for the two columns: 

PvsqKs - Pvsghs = PvsgiKs - hs) < o. (22) 

Here surface pressure pyg of water vapor is set to be saturated, so its 
value is determined by the surface temperature and is the same in 



the two columns considered. It follows from (22) that the first col- 
umn, where condensation took place, cannot remain in static equi- 
librium, as surface air pressure appears to be larger than the column 
weight. The surplus of pressure at the surface causes air to rise and 
all water vapor to ultimately condense as the air rises and cools. 
This illustrates the physical role of water vapor contained in the 
air column as a store of potential energy available for atmospheric 
circulation. 

The global mean store of potential energy per unit atmospheric 
water vapor mass is of the order of Pvs/Pvs ~ RTs/M^ ~ 1.3 x 
10^ J (kg H20)~^ at the global mean surface temperature ^ 
288 K. The global mean precipitation rate is 11 ~ 10^ kg H2O m~^ year" 



(L'vovitch, 1979[ ). Thus, the global rate of potential energy re- 



lease associated with water vapor condensation is of the order of 
URTs/My ~ 4 W m-2. This potential energy flux is sufficient to 
drive the general atmospheric circulation of Earth. The power of 
the latter has been estimated at around ~ 1% of the global power 



of 2.4 X 10^ W m ^ of the absorbed solar radiation (Lorenz, 1967) 



13 



4 Large-scale and compact circulations 



In a stationary circulation condensation of water vapor in the region 
of adiabatic ascent of air masses must be compensated by an inflow 
of water vapor from the hydrosphere to the atmosphere via evap- 
oration. This follows from the general continuity equation written 
in the form of ([S]) or ([9]). If 7 does not change in the horizontal 
direction, the change of 7 in the vertical direction that is due to 
condensation, also turns to zero. 

In a large-scale circulation where pressure gradient forces are 
compensated by friction forces such that horizontal velocity u does 
not change, the maintenance of a constant value of water vapor par- 
tial pressure in the horizontally moving air masses, dp^/dx = 0, is 
achieved via continuous local evaporation. Stipulation dp^/dx = 
describes circulation that is isothermal in the horizontal direction. 
Evaporation flux should be of the same order of magnitude as the 
condensation flux. The latter exceeds the evaporation flux by no 
more than approximately twofold - at the expense of horizontal im- 
port of water vapor evaporated within the region of descent where 
condensation does not take place. The relative value of evaporation 
dictates the direction of air circulation: the air flow in the lower 
atmosphere is directed towards the region of higher evaporation, 
which entails more intense condensation. 

Namely this pattern determines how the biotic pump of atmo- 
spheric moisture functions. Evaporation from the forest canopy of 
natural forests can exceed evaporation from the open oceanic sur- 
face by over twofold. This leads to the appearance of atmospheric 
flow bringing moisture from the ocean to land to compensate for the 



continental moisture loss due to gravitational runoff ( Gorshkov and 



Makarieva 


2006; Makarieva and Gorshkov 2006 


et al. 


2009 




Makarieva and Gorshkov 




2009a 


), Fij 



2007; Makarieva 



. Intense con- 
densation over forest canopy creates a region of low pressure of con- 
tinental scale that sustains the ocean-to-land air flow. The regular 
and persistent atmospheric moisture delivery of the biotic pump 
becomes possible due to the fact that natural forests control the 
processes of evaporation (via transpiration and intercept) and con- 
densation (via production and emission of biogenic condensation 
nuclei). This regular moisture transport does not allow spatial and 
temporal fluctuations of condensation to develop thus preventing 



14 



32DO 
16DO 

SOO 

400 ^ 

200 

IQO 
50 ■ 



■ 32DO I 

■ 16DO 



■ 32DO ■ 

16DO ■ 

8QO ■ 
4DO 

I 

20O t 

IDO ■ 

50 ■ 



(a) 



(b) 



(c) 



lOOO 20DO 3O0D 

X, km 



IDOO 20DO 

X, km 




















































































■ 












































































10 




















































7- 














< 










1 








St 












































■e- 
































t 






































































































































V 










\. 








































1 
















































































k 


f 


fr- 


































































^ 














































































































































i 
































































L 















































Figure 3: Biotic pump of atmospheric moisture in river basins covered by natural 
forests (black symbols) as compared to disturbed and unforested regions (open 
symbols) . Shown is the mean annual precipitation on land associated with large- 
scale regional atmospheric circulation patterns as dependent on distance from 
the ocean. Arrows on the map indicate which regions were considered, 
(a), (b), (c): data for the world's largest river basin covered by natural forests. 
In the basins of world's largest rivers - Amazon and Congo - precipitation is 
twice higher than over the ocean and does not decrease with growing distance 
from the ocean. Precipitation in the basins of the northern rivers that flow from 
the south grow with distance from the ocean proportionally to the increase in 
solar radiation ( Makarieva et al. 2009 ) . In the Ob basin precipitation declines 



in the deforested region near the place where Ob is joined by Irtysh. 
(d), (e), (f): precipitation over non- forested areas decline exponentially with 
distance from the ocean. In panel (f) temperate forests of the North America 
maintain nearly constant precipitation until the deforested region is reached 
(line 4). 
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the formation of hurricanes and tornadoes over the natural forest 
canopy spread over miUions of square kilometers. All these biotic 
pump features disappear in deforested regions, that, via the transi- 
tional monsoon-type regime of intermittent drafts and floods, ulti- 
mately turn to deserts that are deprived of precipitation altogether. 
Fig. 3. 

In a compact circulation, friction forces are negligibly small com- 
pared to pressure gradient forces. This makes it possible for the air 
masses to accelerate and achieve catactrophic velocities observed 
within hurricanes and tornadoes. In this case condensation rate can 
exceed the local evaporation rate as supported by solar radiation 
by several orders of magnitude. The circulation pattern can remain 
stationary if only, after the local store of water vapor is depleted, 
the pressure field moves along to the adjacent area where water va- 
por is still abundant. Having depleted the water vapor in the new 
area, the pressure field moves to the next area and consumes wa- 



ter vapor there, and so on (Gorshkov 1990, 1995; Gorshkov et al. 



2000). Thus, in this case the import of water vapor to the circu- 



lation area occurs dynamically via movement of the pressure field, 
while the consumption of water vapor occurs thermodynamically in 
the course of adiabatic ascent of air masses. 

Hurricane energetics is determined by conservation of the sum of 
kinetic energy of the radial, tangential and vertical air movement 
and potential energy of condensing water vapor on each stream- 
line, as described by the Bernoulli integral. (While Euler's equa- 
tions are non-linear, Bernoulli's integrals that represent the sum of 
kinetic and potential energy, are linear. Therefore, any bunch of 
streamlines (sum of particular streamlines) is also a Bernoulli in- 
tegral.) As friction is negligible, angular momentum is conserved 
along each streamline. The apparent non- conservation of angular 
momentum within the hurricane is caused by a superposition of dif- 
ferent streamlines that correspond to different boundary conditions 
(see Section 8). 

Re-distribution of energy between the radial, tangential and ver- 
tical velocity components of the air flow results in the appearance 
of an eye of hurricane or tornado near the condensation center. Due 
to angular momentum conservation, kinetic energy of the tangen- 
tial flow grows more rapidly towards the center than does kinetic 
energy corresponding to the radial and vertical velocity. Therefore, 
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all energy is ultimately converted to the energy of the tangential 
air flow, while the radial and vertical flow components cease to ex- 
ist (their energy becomes zero). As far as condensation rate and 
the associated potential energy are proportional to vertical veloc- 
ity, condensational potential energy turns to zero as well. All the 
energy store of the Bernoulli integral is now concentrated in the ki- 
netic energy of the tangential air flow. The rotating air flow cannot 
move further towards the center because of the zero radial velocity, 
so a spot of calm weather is formed in the center. In the absence 
of considerable friction, the kinetic energy of the tangential air flow 
can persist for a long time endangering any object that can become 
a source of frictional dissipation. We summarize that the properties 
of hurricanes and tornadoes follow from (1) the properties of the po- 
tential energy associated with condensation, (2) conservation of an- 
gular momentum, (3) conservation of air flow (continuity equation) 
with a high precision, and cannot be understood without taking the 
corresponding physical processes into account. 



5 Condensational pressure gradient in a radially 
symmetrical circulation 

Let us now consider the structure of hurricanes and tornadoes in 
greater quantitative detail. Such compact circulations have a con- 
spicuous condensation center and can be considered in the cylindri- 
cal system. The continuity equation in the cylindrical coordinates 
under the assumption of radial symmetry (no dependence on angle 
(^) have the following form, cf. (|4])-([6]): 

IdNur dNw f dN^ N^dN\ 

+ ^^ = '^\^.. ■ (23) 



r dr dz \ dz N dz ^ 

1 dNdUr ^ dNdW ^ ^ 
r dr dz ' 

1 dN^ur ^ dN,w ^ ^ fdN^ _ N^dN_\ ^^5) 



r dr dz \ dz N dz ^ 

Here u and w are the radial and vertical velocities of the air flow. 



Due to radial symmetry, tangential velocity v does not enter (23)- 



(25). Expanding the derivatives in Eq. (25) and multiplying both 
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parts of the equation by N/N^ we obtain 

1 dur 



dr dz 



^ dw 
r dr dz 



0. 



(26) 



Now expanding the derivatives in the left hand part of Eq. (23 ) and 
using Eq. (26) we obtain the following continuity equation where 



condensation is taken into account: 



dr dr 



dK, 



M^-^^)=^(^-^^)- (27) 



N dz 



Using the equation of state ([T]) we finally have 



u 



dp 
dr 



1 dp^ 
7 dr 



I 1 dp^ X up \ 
-^7(-^--^) (28) 



1 dp 

p dz 



or 



(9 In 7 



(?7 



When the condensation area is horizontally isothermal, we have 
dpy/dr = 0, so Eq. (28) takes the form similar to (11)-(12): 



dp 
dr 



w Ap 

u h-y ' 



Ap=P'J= Py, k 



-1 



k 



-1 



h 



-1 



1 dpy 



h 



-1 



1 dp dp^ 
pdz^ dr 



0, 



(30) 



(31) 



Pv dz 

where scale height hy and h are determined from the conditions of 
adiabatic ascent and hydrostatic equilibrium for moist air including 
the remaining non-condensed water vapor, see Section 3. As in 
(11 )-(12), h.y has the meaning of the scale height where most part of 
water vapor condenses. All variables entering (30) and (31) depend 
on r and z. 

To find the dependence of pressure gradient (30) on r and z 
let us use the continuity equation in the integral form. Air fiow 
converging towards the condensation center penetrates via vertical 
round wall of radius r, circumference 2nr and height and leaves 
the condensation area via horizontal disk of area vrr^: 



2T[rh^uN 



271 / NdWrdr, or 



1 dNur 
r dr 



NdW. 



(32) 
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Here we took into account that all water vapor that entered the 
circumfurence has condensed and it is only the dry air component 
that leaves the condensation area. Recalling the smallness of 7 <^ 1 
and considering that the relative pressure difference between the 
hurricane center and the outskirts does not exceed 10%, one can 
neglect the corresponding changes of ~ A^^ in (32). Then we 
have: 

w 1 dur u du 

(33) 



dr 



u 

r dr 



Putting (33) into (30) we obtain 



dp Ap dur 
dr ur dr 



d\\iur 
Ap— — , Ap 
dr 



PI 



(34) 



where p and 7 are the corresponding values at the surface or, more 
accurately, at the height of boundary layer where relative humidity 
reaches unity and condensation commences; Uc is the condensational 
velocity, which is the velocity scale that determines the magnitude 
of kinetic energy to which potential energy Ap can be converted. 



6 Profiles of pressure and velocities in hurri- 
canes and tornadoes 

Thus, potential energy P{r) that is responsible for the formation 
of hurricanes and tornadoes, depends on radial velocity u{r) (the 
latter itself appears due to conversion of potential energy P(r) to 
kinetic energy) and is equal to 

P{r) = p{r) = Ap In ur + p{rp). (35) 

The Bernoulli integral of Euler equation for a streamline has the 
form 

<l>(r) = -psiu"^ + v'^ + w^) + Aplnur = $(rp), (36) 

where is a fixed radius, v is tangential velocity that is perpendic- 
ular to radius r. As a boundary condition to fix it is natural to 
set u{rp) ~ 5 m s"^ to be equal to the average wind velocity out- 
side the area of the considered compact circulation. Conservation 
of angular momentum constrains the dependence of f on r. Using 
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the continuity equation (33) that relates w and u to each other, we 
then have 



vr = VpTp, V 



Tp W 



1 dur 



r dr 



We now go over to the following dimensionless variables: 



u 



X = 



Ap 



2Ap 2p 



P 



7- 



P 



Ur 



Ur 



w dp 1 dp p 
Uc ' dr Ap dx ' ^ Ap ' 



Ur 



Ap 



a, Vp 



v{rp), Ap 



^p^p 1 -^c 



U^Tp. 



(37) 

(38) 

(39) 
(40) 



Here Uc is the condensational velocity scale (34), Ac is the conden- 



sational angular momentum, a represents the dimensionless angular 
momentum A = Ap (per unit mass density p) that is conserved. 

Leaving the pressure and velocity notations unchanged for the 
dimensionless variables (which is equivalent to choosing the units Uc 
and Ap of velocity and pressure measurements equal to unity) we 
obtain for (36): 



^{x) = + v'^ + w'^ + Inux = $(1), 



a 1 dux 

V = -, w = [5 — — , [5 
X X dx 



p{x) = \nux + p(l). 



(41) 
(42) 



Putting (42) into (41) we finally obtain the following equation on 



radial velocity u: 

a2 + (3^u^ 



u' + 



x^ 



X dx \dx 



+ \nux = $(1). (43) 



Due to the Earth's rotation, in the inertial frame of reference the 
air volumes rotate with an angular velocity u = Qsini), where i) is 
the latitude angle where circulation takes place, = 27r/r, r = 24 h. 
Far from the condensation area due to friction effects the rotation 
is similar to solid body rotation. There is no radial convergence 
of the air flows towards the condensation center. Angular momen- 
tum is proportional to cur^ and declines proportionally to with 
decreasing r. Radial velocity and air convergence towards the con- 



densation center arise at r 



where the pressure gradient forces 
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associated with the on-going condensation significantly exceed the 
friction forces. Starting from r < Vp angular momentum in the iner- 
tial frame proportional to ur^ = VpVp = vr, Vp = uVp, is conserved. 
Tangential velocity increases towards the center in accordance to 



(37) as a consequence of the radial symmetry of the condensation- 
induced pressure gradient forces and the smallness of friction forces 
that can be put approximately equal to zero at r < r^. 

In the system of observations, due to friction between the air 
and the surface, the air is motionless at r > rp and the angular 
momentum is equal to zero. At r < rp there appears a non-zero 
radial velocity and the air masses start to converge towards the 
center. At this moment the Coriolis force that is perpendicular to 
radial velocity vector and to the vector of angular velocity is non- 
central; it starts to curl the converging air masses, which leads to 
the increment of tangential velocity and angular momentum. This 
occurs until tangential velocity becomes considerably larger than 
radial velocity. Then the Coriolis force becomes a central force, and 
the angular momentum is further conserved. 

Tangential velocities in the inertial frame of reference v {v = 
VpVp/r = UTp/r) and in the system of observations Vob are related as 
follows: 

frl \ fl \ 

Vob = V — Lor = uj \ r ] = Vp \ x ] = Vp , x = — . 

\ r / \x J X Tp 

(44) 

Equation (44) describes the fact that in the system of observations 
from the value of tangential velocity in the inertial frame, where 
angular momentum is conserved, one substracts tangential velocity 
equal to the velocity of a solid body rotation ujr. Angular mo- 
mentum per unit mass density p in the system of observations is 
proportional to Vobr: 

r. 



Aob = VobT = wrj ( 1 - ^ ) = A(l - x^), (45) 



where A = urp = VpVp is the angular momentum that is conserved 
in the inertial system. Change of angular momentum in the sys- 
tem of observations as the air moves radially towards the center is 
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determined by Coriolis force: 
dAob 



dt 



■ —2[r[(jj V{r)]] = —2u:ur, V= u + v+ r, (Vr) = ur, 

, , , r dA , r dAdt , f\ 1, 

Aob{r)= / -Trdt= / ——dr = - I 2ujur-dr = 
Jr^ dt Jrp dt dr J^^ u 

dr 

u = Aob{rp) = 0, 



dt 



(46) 



which coincides with Eq. (44). Therefore, at x <^ 1 (r <^ Vp) the 
angular momenta in the system of observations and the inertial sys- 
tem are conserved and coincide. Due to this fact it is convenient to 
consider hurricanes and tornadoes in the inertial system where the 
Coriolis force is absent. Thus, the conserved dimensionless angular 



momentum a is defined as follows, see (40): 



a = aob = a(l — x^) ~ a, x <^ 1. (47) 

Uc 

An important peculiarity of a compact circulation like hurricanes 
and tornadoes is the fact that in order for the circulation to remain 
stationary (independent of time), it is necessary that the entire cir- 
culation pattern moves as a whole with velocity U along the Earth's 
surface. Within each streamline the air becomes deprived of water 
vapor after it ascends in the condensation area. When all the air in 
the condensation area ascends to height h^, the circulation pattern 
has to move over a horizontal distance of the order of the circulation 
size Tp. Consequently, the horizontal velocity of the system move- 
ment is [/ ~ wvp/h^, where w is the mean vertical velocity of the 
ascending air within the circulation pattern. 



In hurricanes the vertical velocity is small, so /3 <^ 1, see (42) 



and can be put equal to zero in (43). Then Eq. (43) assumes a 
simple form: 

,2 



M^ + ^ + ln ( — X ) = a^+M?, p(x)=ln(— x), Ui = u(l). (48) 
x^ \Ui / V^i 



Vertical velocity w is related to radial velocity u by Eq. (42). Point 



of maximum u = Um that is obtained from du/dx = corresponds 
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to X = Xrn = V2a. At smaller x < Xm the value of rapidly declines 
with decreasing x. Contribution of the logarithmic term — \ii{u/ui) 
vanishes as u decreases down to Ui < Um- Further increment of 
— ln('u/'Ui) at M < Ml does not have a physical meaning. Indeed, 
since the water vapor is depleted in the end of the streamline at 
some small x, the condensational potential turns to zero. We thus 
assume that where radial velocity decreases from its max- 

imum value to u = Mo = Ui, is the point where the condensational 
pressure gradient ceases to act (see Section 9). At this point the 
radial and vertical velocities as determined by condensation within 
the considered area are effectively equal to zero. All condensational 
energy is concentrated in the kinetic energy that corresponds to the 
tangential velocity of the air flow. 



7 Rotation of the eye 



According to Eq. (48), the value of Xq at which -u = mq = Mi is a 



root of the following equation. Fig. 4: 



- Inxo = ^(1 - Xg) ^ ^, Xg < 1; xo = xo(a). (49) 

Xq Xq 

Point Xo determines the radius of the hurricane eye for the case 
when the air in the eye were motionless with a zero kinetic energy 
of rotation. However, in this case at x = xo the tangential velocity 
reaches its maximum value vq = a/xQ = uoXo- Therefore, the air 
in the eye cannot be motionless being in contact with the rotating 
windwall at In the absence of friction losses between the 

eye and the windwall, the eye should rotate as well with a constant 
angular velocity Uq determined by the tangential velocity of the 
windwall at x = xq. Energy for the rotation of the eye can be 
only borrowed from the kinetic energy of the rotating windwall, 
because at x = xq the radial and vertical velocities, as well as the 
potential energy of condensation, are close to zero. The energy for 
eye rotation is therefore a certain share of the full potential energy of 
water vapor condensation that was converted to the kinetic energy 
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of the windwal )^r°l 

Transfer of kinetic energy to the eye reduces the maximum tan- 
gential velocity of the hurricane that is observed within the wind- 
wall from fo to Ve = a/x^- The eye radius expands from xq to 
Xe> Xq. The air in the eye rotates with a constant angular velocity 
uie = Ve/xe < ujq and tangential velocity v = UeX, x < Xe- For 
the definitiveness sake, below we continue to call xq as the radius 
of the (motionless) eye while referring to Xe as to the radius of the 
windwall (of the rotating eye). According to the energy conserva- 
tion law, the kinetic energy of eye rotation should be equal to the 
kinetic energy of rotation of the windwall segment Xq < x < Xe with 
tangential velocity v = a/x which would take place if the eye re- 
mained motionless, see Fig. 5a,b. Neglecting the air density change 
in the eye as compared to the windwall (see footnote 5) and can- 
celling the common multiplier 27r, this stipulation can be written in 

^The pressure gradient force is perpendicular to isobars. Coriolis force is perpendicular 
to the velocity vector and is proportional in magnitude to velocity. If the two forces act in 
the opposite directions and compensate each other, one can determine the constant value of 
velocity corresponding to movement along the isobars. In the inertial frame of reference for 
the case of radially symmetrical closed isobars this corresponds to the equality between the 
centrifugal force that depends on constant tangential velocity v and the pressure gradient force: 
v'^/r + dp/dr = (the so-called cyclostrophic balance). If there is no friction, the angular 
momentum is conserved, and the rotational movement with zero radial velocity and constant 
tangential velocity can continue infinitely, similarly to satellite rotation on a terrestrial orbit. 
Friction decreases tangential velocity and, hence, the centrifugal force; this disturbes the 
balance of forces and causes some radial movement towards the center of the isobars. In the 
existing theories of hurricane formation the pressure gradient force is approximated from the 
condition of approximate cyclostrophic balance using the dependence of tangential velocity 
on radius and the fact that the observed radial and vertical velocities are small compared to 
the former. Friction is con sidered as the major cause of the convergence of air masses towards 
the hurricane center (e.g., [Smith et aL||2008"] l. The cause-and-effect link between the radial 
convergence and the existence of the radial pressure gradient is neglected. However, nature 
works the other way round: it is namely the radial convergence of moist air with a non-zero 
radial velocity u and the physically inseparable ascent of air masses with vertical velocity w, 
see \33\ , that determine the pressure fall from the outer environment towards the hurricane 
center, see ( |34[ l. If there is no radial convergence and the radial velocity is zero, then there 
is no ascent of moist air, no condensation of water vapor, no condensation-induced pressure 
gradient, no close isobars and no approximate cyclostrophic balance - in short, no circulation 
and, consequently, no friction. This is the major difference of the condensational potential 
from the velocity-independent gravitational potential that keeps the rotating satellite on its 
orbit. The condensational potential energy as the air converges towards the center is mainly 
converted to the kinetic energy of rotation determined by tangential velocity v. This results 
in the development of huge velocities namely due to the absence of friction. Friction can only 
lead to dissipation of this energy impeding the air acceleration and preventing the development 
of maximum possible velocities. 

^'^The appearance of eye rotation due to the contact of the eye with the rotating windwall 
represents a re-distribution of the ordered kinetic energy. It cannot not be considered as 
turbulent friction losses (i.e., energy losses on the formation of the chaotic turbulent air flows). 
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Figure 4: The regime of hurricane existence corresponds to the two solution 
branches of equation ( 49 1 . The green curves in all panels correspond to the 
branch where the eye radius xo(a) grows up to a;o — \Jla — 0.6 as a increases 
up to a « 0.43. The red curves in all panels correspond to the branch where the 
eye radius grows from xq = 0.6 to Xq = 1 with a diminishing from 0.43 to zero. 

(a) : Two branches of the solution of the approximate equation ( 49 1 . 

(b) : Maximum dimensionless tangential velocity i'e(a) 



(51), 



ajxe = a/(1.28a;o) 
The branch with lower velocity corresponds to widening of the eye accom- 



panied by the decrease of condensational energy. 

(c): Absolute values of maximum tangential velocities Ve.(a) = u^ajxe depend- 
ing on temperature for five different values of parameter TpSini?. For example, 
Tp sin = 300 km corresponds to a streamline with the outer radius — 600 km 
formed at latitude -d — 30°. The value of Uc grows exponentially with increas- 
ing temperature as = [27(r)i?r/M]i/^ see ([ss]), 7(r) = Pv{T)/p, p^{T) 
conforms to the Clausius-Clapeyron equation dp^/py = {L/RT)dT/T, sec (16), 

7 = 0.042 at r = 303 K. 

Black dots represent the data from Fig. 3 of Michaels et al. 



(2006) that corre- 



spond to maximum wind velocities observed in 270 tropical cyclones. 
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the following form, see Fig. 5a: 



X xdx 



xo 



xdx 



x^ 



xl 



(50) 



Performing the integration and cancelling the common multiplier 
we obtain 



1 

4' Xo 



.0.25 



1.28; Xe{a) = 1.28a;o(a), Ve{a) 



In 

(51) 

Tangential velocity v grows with decreasing x aX x > Xe and declines 
with decreasing x at x < Xg as follows, see (47): 

a 



V 



< 



X 

ax 
xl ' 



X > Xp, ux 



X < Xe, fX 



Mr, 



const., Lo = Vt sin'i?, 



(52) 



ax 
x^ 



7^ const. 



The dependence of the maximum value of tangential velocity Ve = 
v{xe) on a is shown in Fig. 4b. 

Stationary existence of the rotating windwall and the rotating 
eye modifies pressure profile p{x) at x < Xg in such a manner that 
the pressure gradient force and the centrifugal force within the eye 
coincide. In order that an additonal pressure drop could form within 
the eye where no condensation takes place, some part of the air must 
be exported away from the eye outside the hurricane area (x > 1) 
during the eye formation in such a manner that the parts of the 
pressure profile inside and outside the rotating eye joined smoothly 
and featured no discontinuity at the windwall x = Xg. This condition 



can be written as (Makarieva and Gorshkov, 2009b) 

,2 



dp 
dx 



V 2 

2 — = 20; X, p[x) 

X 



+p(0). 



Up 



a 
xl 



p{x) = In 



u 



-X 



Ui 



, p{l) = 0, X > Xe 



X < Xe; (53) 



(54) 



(Multiplier "2" in (53) arises from the equation dp/dr = v'^/r as one 
goes to the dimensionless variables (39): dividing both parts of the 
equality by Ap = psul/2 and changing to the dimensionless velocity 
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v/uc-) In Eq. (54) the condensational potential is set equal to zero 
outside the condensation > 1 (r > Tp). Thus, p{x) describes 

how air pressure declines from its normal value it has outside the 
hurricane. Accordingly, — p(0) = 5p is equal to the total (dimen- 
sionless) pressure fall within the hurricane including the pressure 
fall associated with the eye rotation. In ordinary units, the pressure 
fall is equal to 6pAp = jpSp. Combining (53) and (54) we have, see 
Fig. 5c, d: 



p{x) 



In 



-X 



u 

Ui 
2 2 



+ In I — X 

.^1 



u 



X> Xf. 



X < Xp 



(55) 



Sp 



In 



Up 



-X, 



Ui 



(56) 



The first and second terms in (56) represent the pressure fall within 
the eye and outside the eye, respectively. At a = 0.12 (r^ = 400 km. 



d = 20°, Mc = 83 m s"^ at 7 



Xe ~ 0.15, see (51); Up/ui 
we obtain Sp ~ 3, Fig. 5d 
SpAp/p ~ 37 ~ 0.13 at 7 ^ 



i 0.04) we have xq ^ 0.12, see Fig. 4a; 
2, see Fig. 5b. Thus, from Eq. (56) 



For the absolute pressure fall we have 
0.04, which means that the air pressure 
in the eye center is at maximum 13% lower than the air pressure 
outside the hurricane. 



Note that the two branches of pressure p{x) (55) inside and out- 
side the eye join smoothly at x = Xe, while the derivative dp/dx 
features a minor discontinuity at this point. This discontinuity is 
practically unnoticeable in Fig. 5d. This discontinuity has a clear 
physical meanining, but it is difficult to observe it at small a. If 
using the smallness of a, Xg and x^ one chooses Xg = Xm at the 
point where radial velocity has its maximum, then the derivatives 



dp/dx also coincide in this point, see (49). This is a consequence of 
Eq. ( 48 ) and of the fact that du/ dx = Q a,t x = Xm ( Makarieva and 



Gorshkov, 2009b) 



function of a (|51[), xq = 
0.43 



Equation (49) represents an equation on the value of xq as a 

xo(a). Fig. 4a. The value xq = x^ = 0.6 
0.36 and the approximation 



is reached at a 



where Xq 



(49) still holds. Further widening of the eye radius xq involves the 
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Figure 5: Dimensionless radial u (48) and tangential v (52] velocities and pres- 
sure fall p{x) (55) as dependent on radius to the hurricane center x. 
(a), (b): velocities and pressure assuming the eye is motionless; 
(c), (d): velocities and pressure with an account of eye rotation; 
Xq is eye radius (the eye is motionless); 
Xe is windwall radius (the eye is rotating); 
Xm is the point where radial velocity reaches its maximum; 
a = TpLo/uc is the dimensionless angular momentum that is conserved on the 
streamline outside the eye; Vp is the outer radius of the hurricane streamline; 
Uc — {2/S.p/ pY^'^ is the condensational velocity scale, see (38) and notes to 
Fig. 4c; Ap = is pressure drop due to condensation, 7 = Pv/p is the relative 
share of water vapor, p is air pressure outside the hurricane. 
The curves correspond to the following parameters: — 400 km, u — rtsin-d = 



2.5x10"'^ s 
a = 0.12, Uc 
ui = 0.06. 



, d = 20°, 7 = 0.042, T = 30° C, p = 1.22 kg m'^ p = 10^ J m'^ 
83 m s~^. The black horizontal line in panels (a) and (c) denotes 



28 



second branch of the solution of Eq. ( 49 ) shown in Fig. 4a. For this 
branch the growth of Xq to Xq = 1 corresponds to decreasing a. On 
this branch potential energy arising from condensation that occurs 
in the area Xq < x < 1 and disappears at Xq = 1, a = 0, should 
have been spent on rotation of a very wide eye at low tangential 
velocities, i.e., the hurricane cannot form. 

Transition from this branch with Xo > Xm to the branch xq < Xm 
corresponds to formation of a hurricane. In Fig. 4c the areas of 
possible maximum tangential velocities are shown for several val- 
ues of Tp sin as dependent on oceanic surface temperature T, that 
dictates the value of condensational velocity Uc (38). We note two 
peculiarities: (1) temperature that corresponds to the transition be- 
tween the two regimes decreases with diminishing rpsmd and (2) 
for any given temperature, the difference between velocities of the 
two branches grows with decreasing Tpsmd. In the equatorial zone 
where smd — )• 0, at observed surface temperatures the difference 
between velocities on the two branches tends to infinity. Hurricanes 
cannot form. 

As one can see from Fig. 4c, for any given temperature, the 
distance between the branches increases (while the probability of 
jumping from one branch to another consequently decreases) with 
decreasing radius and/or latitude. For example, for curve 5 at 300 K 
the distance between the green and red branches is nearly 80 m s~^, 
while for curve 2 it is about two times smaller. The transition point 
between the two branches for curve 5 corresponds to T = — 5°C, 
which is a temperature when the hurricanes do not form. The tran- 
sition point moves to the region of lower temperatures with growing 
TpSmd. On the other hand, with increasing radius and/or latitude 
the area of allowed velocities becomes narrower. For example, for 
curve 1 at the observed temperatures only a narrow range of veloc- 
ities in the vicinity of 46 m s~^ can be developed. These physical 
limitations shape the observed distributions of hurricane wind ve- 
locities over temperature. Fig. 4c. 

Finally, when the angular momentum a decreases (due to either 
decreasing radius Vp of the condensation well as due to de- 

creasing angular velocity u at low latitudes), the eye radius decreases 
as well. The energy ot the eye rotation, which in the dimensionless 
variables (39) is equal to a^/4, see (50), becomes small. In the result, 
there appears a possibihty for tangential velocity in the windwall to 
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grow to catastrophic values that are characteristic of tornadoes. 



8 Angular momentum: the effect of superposi- 
tion of several streamlines 



Consideration of the energy budget of every streamhne starts from 
r — Tp and ends with r = 0, with the account made for the sub- 
straction of energy to sustain the eye rotation. Eye radius x^. for 
each streamhne depends on angular momentum a = Tpu/uc, i.e., it 
is also determined by the value of r^. Angular momentum of the 
rotating air volumes in the eye changes with distance according to 
ar'^/r^ _ ax'^ /xl- Thefore, the complete information about a given 
streamline is contained in the value of angular momentum a that is 
conserved in the region outside the eye x > x^- 

Hurricanes apparently comprise a large number of individual 
streamlines, some of which can be easily traced on the satellite im- 
ages as individual "tails" spiraling towards the hurricane center - 
the eye. These streamlines can start at various radii and possess 
various values of angular momentum a that is conserved outside the 
eye. 

The average angular momentum of the hurricane can be obtained 
considering the sum of kinetic energies of rotation for all streamlines 
as dependent on distance r to the circulation center. For each i-th 
streamline starting at r = Vpi the dependence of angular momentum 
ai{x) and tangential velocity Vi{x) on relative distance x has the 
form 



ai{x) = ai 



X 



■d{Xei -X)+ d{x - Xei)'d{Xi - x)) 



, Vi{x) 



OiiX) 
X 



Xi 



' pi 



' pmax 



X 



< 1, Oj 



' pmax 



Ur 



1, ifX>0 
0, if X < ■ 
(57) 



Here x is defined relative to the streamline with maximum ultimate 

radius Vp^ax ^ ^pi and maximum conserved angular momentum; 
'd{X) is the step function. Mean angular momentum a{x) as depen- 
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dent on x can be obtained from the following relationship: 



where aj(x) and fj(x) are defined in ([57|). The last approximate 



equality in (58) is written from the condition ^jaf(a^) ~ no^{x) 
a^{x) ~ a^^xjoi^x) to qualitatively evaluate the behavior of a{x). As 



one can see from Eq. (58 ), the values that enter the sum diminish 
with decreasing x. Hurricane eye acquires a complex structure that 
is formed by different streamlines at different radii Xei. The closer to 
the center, the smaller values of angular momentum are that form 
the hurricane eye. These two facts testifiy that, in agreement with 
observations, a{x) decreases towards the center despite the angular 
momentum is conserved within each particular streamline. 



9 The pole approximation 



The circulation system of hurricanes and tornadoes can be consid- 
ered in the pole approximation. Such a consideration preserves the 
main physical features and the quantitative estimates of velocities 
and condensational potential. In the pole approximation one pre- 
serves the terms that display most singularity at a; 0, namely, the 
major powers 1/x and Inx. This is equivalent to discarding du/dx 



as being small compared to u/x in the Bernoulli integral (43) and in 



the Euler equation for the streamline that corresponds to differen- 



tiating both parts of (43) over x. In this case the Bernoulli integral 



(43) simplifies to become (Makarieva and Gorshkov, 2009b) 



+ v'^ + w'^ + Inx = 0, 



2 02"^^ 
— , w = /3 ^, 
x-^ X'^ 



u 



Inx 



X 



, p{x) = Inx. 



(59) 
(60) 



For hurricanes, where /3 ^ 1, Eq. (60) takes the form 



u 



\nx 



x^ 



(61) 
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The point of maximum Xm = y/2a. For maximum radial velocity 
uL we have 



\ii{V2a) - l- (62) 



Velocities u and w become strictly equal to zero at x = xq defined 
from the equation 

— Inxo 2 ~ (^'^) 

Xq 

Tangential velocity reaches its 

As one can see from the comparison of (41) and (59), the pole ap- 
proximation consists in replacing condensational potential ln[(M/Mi)x 
for Inx, i.e., it does not depend on the initial velocity U\. The gra- 
dient of the potential 1/x -|- vT^dujdx is replaced by 1/x, i.e. the 
term containing du/dx is discarded. This changes the profile of the 
potential in the region from x = 1, m = Mi = u{l) to x = Xq, 
u = uq = m(xo), while the change of the potential (pressure dif- 
ference 6p) is not affected by the approximation. As far as with 
decreasing x the value of u first grows at x > x^ = V^a and then 
declines at x < x^, the exact potential declines more slowly than 
the potential in the pole approximation at x > x^ and more rapidly 
than the latter at x < x^- 

Tornadoes, where /3 > 1, preserve the main features of the pole 
approximation. Angular momentum in the tornado is very small, 
a = TpUj/uc > h^uj/uc < 0.004 (at ~ 10 km, a; ~ 3 x 10~^ s, Uc ~ 
80 m s~^). Due to the fact that the eye in tornadoes is extremely 
narrow, Tq = h^Xo = ~ 20 m, the eye rotation practically does 

not take any energy from the streamline. Neglecting the difference 
between xq and Xe we have from Fig. 4a that Xq ~ Xq ~ a/2, see 
Fig. 6a. From (60) we have 6p ~ p{l) —p{xo) = — Inxo ~ 6, Fig. 6b, 
so that 6pAp/p ~ 67 ~ 0.3. In the center of the tornado the air 
pressure can drop by 30% compared to the outer environment. This 
causes tangential velocity to grow to very high values of the order of 
{a/xo)uc ~ 2uc ~ 160 m s~^. Fig. 6a. In contrast to the hurricane, 
while approaching x = Xq the radial velocity u in tornado decreases 
proportionally to x^, see (60) and Fig. 6a, and has a maximum at 
X = 0.69 (—Inx = 1/2.). Vertical velocity w (59) in tornado has the 
same behaviour as radial velocity u in the hurricane, see (60). In 
the point of maximum Xm = V^a the vertical velocity reaches Wm = 
Uc{—ln(\/2a) — 1/2)^/^ ~ 160 m s~^, see (62). This maximum value 
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Figure 6: Wind velocities and pressure profile in a tornado calculated in the 
pole approximation. 

(a) : vertical tangential v and radial u velocities in units of Uc = {2Ap/ p)^/"^ 
( 38 ) depending on distance x to the condensation center (x measured in units 
of the condensational scale height (11) at /3 = 1) with dimensionless angular 
momentum a = h^co/uc, xq is the eye radius that satisfies Ihxq = Ix\, Typical 
values at 30°C are 80 m s-\ - 10 km, Fig. 2, Ap - 4 x 10^ Pa, 
o = 0.004, xo = 0.0016. 

(b) : pressure fall p(a;) in units of Ap from the outer edge of the condensation 
area at a; = 1 to the windwall - eye radius at x = xq. 



of vertical velocity coincides with the value of tangential velocity 
reached at xq ~ a/2, Fig. 6a. Radial velocity u in tornadoes changes 
much more slowly at small x than it does in hurricanes, du/dx <^ 
u/x, which makes the pole approximation more exact for tornadoes 
than it is for the hurricanes. 

We note one more peculiarity of the tornado. The eye radius is 
xo ~ a/2. Fig. 4, a = h^u/u^, = (2Ap/p)i/2 = (2p/p) 1/2^1/2^ gee 
(34). The water vapor relative partial pressure 7 decreases approx- 
imately exponentially with growing height and practically turns to 
zero at 5-8 km, see (18). Therefore, the eye radius in the tornado 
grows proportionally to 7"^/^ with increasing height z. This creates 
a conspicuous "mouth" of the tornado funnel - the funnel is narrow 
at the surface and widens as the height grows. This is another indi- 
cation that tornado represents a single streamline in contrast to the 
hurricane, where the eye is formed as the result of many streamlines 
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being superimposed. 



10 Major features of atmospheric circulation in- 
duced by water vapor condensation 

1. The physical cause of the condensation-induced circulation con- 
sists in the fact that hydrostatic equilibrium of the moist terrestrial 

atmosphere is unstable and cannot exist (Sections 1, 2). 

In an isothermal atmosphere partial pressure of water vapor would 
exponentially decline with height decreasing twofold per each nine 
kilometers of height increment. According to Clausius-Clapeyron 
law, saturated water vapor would follow the same exponential dis- 
tribution if the air temperature dropped by approximately ten de- 
grees K per each nine kilometers of height increment, i.e., with a 
lapse rate of 1.2 K km~^. However, adiabatic ascent of a volume of 
moist air causes air temperature to drop at a significantly higher rate 
bounded between approximately 4 and 9.8 K km~^. Consequently, 
a random adiabatic displacement of moist air upwards makes the 
excessive water vapor condense and leave the gas phase. Pressure 
in the adiabatically ascending air drops. There appears an upward- 
direction evaporative-condensational force that sustains continuous 
adiabatic ascent of the moist air and causes air to circulate in the 
vertical and horizontal directions. 

2. Stationarity of the annually averaged state of the atmosphere 
and the surface testify to the equality between the global mean 
rates of evaporation and condensation rates. However, evaporation 
is maintained by the stationary flux of solar radiation. Mean power 
of the evaporation flux cannot exceed the power of solar radiation. 
In contrast, the power of the flux of water vapor condensation is 
not related to the solar radiation power and can locally exceed the 
evaporation power by one-two orders of magnitude (e.g., in hurri- 
canes and tornadoes) or become by the same factor less intense than 
local evaporation (in the regions of the anticyclonic descent of air 
masses). These contrasting properties of evaporation and conden- 
sation are crucial for the condensation-induced circulation. 

3. In a large-scale horizontal circulation the air flows from the 
donor area to the acceptor area in the lower atmosphere and in the 
opposite direction in the upper atmosphere. In the donor region the 
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air masses descend; there is no condensation, but the evaporation 
is present. In the acceptor area the air masses ascend; there is 
condensation; both locally evaporated water vapor and the water 
vapor imported from the donor area undergo condensation. The 
rate of condensation is thus always higher in the acceptor area than 
in the donor area; this sustains the pressure gradient necessary for 
the maintenance of the air flow. Water returns to the donor area 
in the liquid phase. In a large-scale circulation total evaporation 
can coincide with total condensation (then the circulaton can be 
stationary in the long-term). However, the circulation itself is only 
possible due to the absence of equality between local evaporation 
and condensation (they do not coincide in either donor or acceptor 
area). Condensation rate in the acceptor area must be higher than 
in the donor area. This is stably achieved by higher evaporation 
rate in the acceptor compared to the donor 

4. In large-scale (Hadley, biotic pump of continental-scale natu- 
ral forest) and medium-scale (moonsoon) circulation systems friction 
forces prevent significant acceleration of air masses and do not allow 
high wind speeds to develop. If there is no control over condensation 
and evaporation, the areas where such circulation systems develop 
are also prone to the appearance of frequent floods, droughts, torna- 
does and hurricanes of varying intensity. The condensation-induced 
circulation can be regulated only with use of the genetically based bi- 
ological programs of natural forest ecosystems, where the necessary 
soil moisture store is continuously maintained by biotic regulation 
of evaporation and condensation. 

5. Compact circulation takes place where friction is small. Con- 
densation power that locally exceeds the power of evaporation by 
many times accelerate the air to catastrophic velocities observed 
in hurricanes and tornadoes. In a compact circulation there is a 
pronounced condensation center, to which the rotating air masses 
converge. The condensation center moves continuously in the hori- 
zontal plane. This movement guarantees that the circulation system 
is supplied by new amounts of water vapor which is depleted as the 
system functions. If compact circulation do not move, they dissipate 
(similar to the movement of an animal over its feeding territory). 

6. With increasing mixing ratio 7 of water vapor in the moist 
air and as this ratio approaches unity due to the removal of the 
non-condensable air constituents (oxygen and nitrogen), the moist 
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adiabatic lapse rate approaches the hydrostatic value of 1.2 K k,~^. 
The evaporative-condensational force that supports the adiabatic 
ascent tends to zero. The adiabatic ascent stops; the atmosphere 
becomes vertically isothermal; condensation-induced circulation sys- 
tems of all types disappear. But with increasing temperature the 
evaporative-condensational force grows exponentially. 

7. The condensation-induced circulation arises on horizontally 
isothermal surface and does not demand an external differential 
heating to arise. 
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